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Abstract 

We investigate the case of [Cll] 158/zm observations 
for SPIC A/SAFARI using a three-dimensional magneto- 
hydrodynamical (MHD) simulation of the diffuse interstel- 
lar medium (ISM) and the Meudon PDR code. The MHD 
simulation consists of two converging flows of warm gas 
(10 4 K) within a cubic box 50 pc in length. The interplay 
of thermal instability magnetic field and self-gravity leads 
to the formation of cold, dense clumps within a warm, 
turbulent interclump medium. We sample several clumps 
along a line of sight through the simulated cube and use 
them as input density profiles in the Meudon PDR code. 
This allows us to derive intensity predictions for the [Cll] 
158^m line and provide time estimates for the mapping of 
a given sky area. 

Key words: Galaxies: formation - Stars: formation - Mis- 
sions: SPICA - macros: M^X 



2. The MHD simulation 

To model the turbulent interstellar medium, we used a 

data cube from an MHD simulation performed by Hennebelle et al 



1. Introduction 

The fine structure line of ionised carbon [Cn] at 157.7 /im 
is one of the most prominent far infrared (FIR) line. It is 
widespread in the Milky Way as revealed by the C OBE- 



FIRAS observations, and detected in the diffuse ISM (Ingalls 
12002 ). photo- dissociation regions a s well as local and dis- 
tant galaxies (jMaiolino et all 12009 . and references therein). 
[Cn] is one of the main cooling lines of the neutral ISM, 
in all environments where carbon is mostly ionised, that is 
where UV-photons with hv > 11.3 eV are available. Unfor- 
tunately, the present information on the spatial structure 
of the [Cn] emission is very limited, because of the small 
telescope size of previous far infrared missions. While Hcr- 
schel and SOFIA will improve the situation for bright re- 
gions, their sensitivity will be limited for the extended 
emission from either the diffuse interstellar medium, or 
the outskirts of molecular clouds. The structure of the 
[Cn] emission is expected to reveal the underlying struc- 
ture of the matter, and provide information on the mecha- 
nisms responsible for the fragmentation. In this paper, we 
explore the perspectives offered by SPICA for mapping 
the interstellar medium. 



(2008|) with the RAMSES code (|Tevssierl . l2002HFromang et all 
120061 ) , which makes use of adaptive mesh refinement (AMR) 
methods. The cube is 50 parsecs on each side, with poten- 
tially 14 levels of refinement, leading to a maximal resolu- 
tion of 3 mpc or 630 AU. However, to simplify the analysis, 
we used a regularly gridded (1024 3 ) version of the cube, 
with a pixel size ~0.05 pc. 

In the initial setup of the simulation, two converging 
flows of warm gas (8000 K) are entering the cube along 
the X axis from two opposite sides. The flows' velocity 
is twice the sound speed of the warm neutral medium 
(WNM), with transverse and longitudinal modulations of 
100% amplitude with a lOpc periodicity. Periodic bound- 
ary conditions are imposed on the remaining four sides. 
The simulation includes a consistent treatment of hydro- 
dynamics, magnetic fields, atomic cooling processes and 
self-gravity, and therefore implies some heavy computing 
(-30,000 CPU hours). 

FigureQ]shows the total gas column density in the cube 
when viewed along the X direction, which we take as the 
observer's direction to ensure statistical homogeneity of 
the observed field. Our chosen line of sight is displayed, 
and passes through one of the cold dense clumps formed 
by thermal instability in the turbulent flow. These clumps 
6l ai.10 3 cm' 3 and T — 10 K) are found along filamentary 



structures amidst a much more diffuse and hot interclump 
medium (n ~ 1 cm -3 and T — 10 4 K). 

Figure[5]shows a total gas density cut through the cube 
in the X — Y plane containing the line of sight. Complex 
density structures along this line appear quite clearly. 

When it comes to simulating the observation of this 
field, we shall need to specify the angular pixel size. To 
match the SAFARI specifications, let us consider that the 
cloud is located 1.75 kpc away. It is therefore 1.6° across, 
and the pixel size is 5.75", i.e. that of the SAFARI pixels. 

3. PDR code 



The Meudon PDR code (|Le Bourlotet al.l . ll993ULe Petit et al 



20061 ) is a publicly available modeQ aimed at describing 



1 http://pdr.obspm.fr 
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Figure 1. Total gas column density along the X direction 
- that of the converging inflows - in the MHD simulation. 
The figure on the right shows the entire 50 pc x 50 pc 
field. The figure on the left presents a zoom on a clump at 
the top right corner of the field. The line of sight used in 
the analysis is marked with a black cross. 




Figure 2. X — Y total gas density cut through the clump of 
interest. The line of sight used in the analysis is indicated, 
and the observer is located at the right side of the figure. 



the atomic and molecular structure of interstellar clouds, 
especially photon-dominated regions (PDR). It considers 
a plane-parallel slab of gas and dust illuminated on either 
or both sides by an incoming radiation field which can be 
the interstellar standard radiation field (ISRF) or the light 
from a nearby star. The density profile within the slab can 
either be constant, computed from an isobaric constraint 
or user-defined. Radiative transfer in the UV is solved at 



each point in the cloud (jGoicoechea fc Le Bourlotl . 120071 ). 
taking into account continuum absorption by dust grains 
and molecular absorption lines. Thermal balance is com- 
puted via the inclusion of a number of heating (photoelec- 
tric effect, cosmic rays) and cooling processes (infrared 
and millimeter emission lines) and an extensive chemical 
network. With abundances and level populations of the 
most prominent species available, line intensities and col- 
umn densities can be obtained. 

Convergence issues make it difficult to run the code re- 
liably on the more diffuse regions. Consequently, we opted 
not to use the entire line of sight as input density profile 
in the PDR code, but rather extracted dense "clumps" 
(see Fig. [3]) using a threshold value of 50 cm -3 - which 
is about the critical density for C + in electronic collisions 
at T e = 10 4 K - and ran the PDR code on these clumps 
separately. 




Figure 3. Total gas density profile along the line of sight. 
Regions where density exceeds 50 particles per cubic cen- 
timeter are designated as "clumps" and numbered from 
1 to 3. The whole region marked '1-3' is designated as 
"cloud" in the text. 



The one-dimensional geometry imposed by the code 
makes shielding an important point to remember when 
considering results. The complexity of density structures 
in the simulation cube implies that the radiation field 
should be higher in the interior regions than expected from 
the code, as it can be expected that the photons penetrate 
deep into the clouds, and, consequently, that the addition 
of results from separate clumps might give a better idea 
of the output emission. 

To estimate this effect, we also ran the PDR code on 
the entire "cloud" marked T-3' on Fig. [31 so in total four 
times. The density in the small interclump regions does 
not fall to such low values that the code may not converge. 



4. Results 

As explained before, we ran the PDR code a total of four 
times : once for every clump 1 to 3, and once over the 
entire cloud 1-3 (marked in red on Fig. [3]). Results are 
summarized on Fig. [4] and in Tab. [1] 

The effect of shielding due to the one-dimensional ge- 
ometry assumed by the PDR code appears quite clearly, 
as the [Cll] emission from the inner regions of the cloud 
is almost completely suppressed compared to results ob- 
tained when running the code on separate clumps. This 
is of course because photons that are energetic enough 
(hv > 11.3 eV) to ionize atomic carbon are absorbed on 
the surface of the cloud and cannot reach these inner re- 
gions in the one-dimensional approach. 

More realistically, given the complex structures in three- 
dimensions, UV photons may penetrate the cloud much 
more deeply from other directions, and as a result, the 
radiation field should be higher than expected from the 
results of the PDR code applied to the entire cloud. 
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Figure 4- Local emissivity of the [Cll] 158/j,m line obtained 
from the Meudon PDR code, along the line of sight. Each 
plot shows the emissivity computed when using the corre- 
sponding clumps from Fig. as input density profiles to 
the PDR code. 



The line being optically thin (jNakagawa et al. , 1995, 
see also optical depths in Tab.[T]), we can run the code on 
clumps and add up contributions from these to obtain a 
more accurate estimate of column densities and line emis- 
sivities. 

Table 1. Hi column densities, [Cll] 158 fiin integrated 
emissivities and line-center optical depths computed by the 
PDR code for clumps 1, 2, 3, as well as for clump 1-3 
(marked in red on Fig. [3J). The sums of column densi- 
ties, integrated emissivities and optical depths for all three 
clumps are also indicated. 



Clump # 


Nhi (cm" 2 ) 


I c+ (erg.cm 2 .s 1 .sr 1 ) 




1 


7.87 10 19 


6.55 10 -6 


0.39 


2 


1.39 10 20 


5.41 10~ 6 


0.34 


3 


1.57 10 20 


6.86 10" 6 


0.99 


1-3 


1.70 10 20 


7.21 10~ 6 


1.19 


1+2+3 


3.70 10 20 


1.88 10" 5 


1.72 



The data in Tab. [T] hints at an empirical linear corre- 
lation between Hi column density and [Cn] 158/im inte- 
grated emissivity from the cold and dense regions of the 
cube, which reads 



10 6 ergs. cm 2 .s 



5 x 



(Nhi) 



.ST 



Obviously, this is a very crude approach, as one of the 
five data points is actually derived from three others, and 



in any case this relationship may well not apply to more 
diffuse regions. To m ake a compari s on wi th observational 
data, we may refer to llngalls et al.l (|2002f ). who give 



I cm 



ergs, cm 



0.32 x 
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2.58/ 6 



MJy.sr 



using ISO observations of high Galactic latitude clouds 
and IRAS far infrared data. The FIR-Hi correlation is in 



turn given by iBoulanger fc Perault.1 (|1988| ) , using IRAS 
data and combined Hi surveys. Averaged over mid-range 
Galactic latitudes 27.5° < \b\ < 32.5°, we have 



'60 



2.3 10 



-21 



MJy.si— 1 
By combining the two we get 

{c+ 

10~ 6 ergs.cm _2 .s _1 .sr _1 



and 



= 0.18 



'100 



= 0.49 x 



N. 



HI 



10 20 cm" 



which is a factor ten lo wer than our re l ation , and in agree- 
ment with COBE data (jBennett et al.l . ll994l ). On the other 
hand, the [Cll]/Hi ratio derived from BIRT data on the 
Galactic Hn region W43 yields 



l c+ 



10 -6 ergs.cm~ 2 .s 



37 x 



Nhi 



This shows that the [Cii]/Hi ratio strongly depends on 
the physical conditions, and that more work is definitely 
needed to understand this correlation. For the sake of sim- 
plicity, and given that our relationship happens to be a 
compromise between the observational data above, let us 
assume that we can use it to straightforwardly derive a 
[Cll] emission map over the entire field. 

It remains that what we have in the data cube is the 
total gas content, so we need to estimate the atomic frac- 
tion , which we do crude ly using a fi10 to the data in Fig. 
6 of iSavage et ail (|l977h . 

This allows us to derive a histogram of [Cn] 158 /j,m 
line intensities emerging from the cube, using the solid 
angle value dft=6.1 x lO" 10 sr for the SAFARI pixels. This 
histogram is showed on Fig. and we will use it in the 
next section to estimate a mapping speed for the entire 
50 pc x 50pc cloud. 

Before that, however, we should estimate the [Cn] emis- 
sion associated with the diffuse regions of the warm ionized 
medium (WIM) phase, if only to check that this contribu- 
tion can be neglected. To thi s end, we used the cosecant 
law given bv lRevnoldd (|l992lh 



J™™ = 3 6 x i(r 7 csc|&| ergs.cm^.s^.sr- 



To be exact, we use model 2 from iBlackl (|l975l ). although 
the mean total gas density in our cube is (n to t) — 5 cm -3 , 
which would point to model 1. The difference between the two, 
as far as the derived [Cll] intensities are concerned, is negligi- 
ble. 
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Figure 5. Histogram of the [Cll] 158(im fluxes on the entire 
50 pc x 50 pc field, using our empirical relationship and 
the total gas to Hi conversion described in the text. 

Given t hat there is another cosecan t law for the Hi column 
density ( Dickey fe Lockmanl . Il990h . 



(N HI ) = (3.84 esc |6| - 2.11) x 10 20 cm" 2 , 

we may combine the two to have an estimate of the WIM's 
contribution to [Cn] emission, 



(J WIM } 



io- 



ergs.cm 



-2 



.sr 



0.094 



HI 



10 20 cm" 



2.11 



Using the same total gas to Hi conversion as before, this 
time on the mean column densities over the whole field, 
we find that (N to t) = 8.1 10 20 cm" 2 implies a molecular 
fraction of about 15%, so that (N HI ) = 6.8 10 20 cm" 2 
and therefore (/$+ M ) = 8.3 x 10" 7 ergs. cm 2 .s 1 .sr 1 , or 
(X^ M ) = 5.1 x 10" 19 W.m" 2 in terms of emerging flux 
for each pixel in the map. Consequently, this contribution 
to the total [Cll] emission is completely negligible. 



5. SAFARI MAPPING SPEED 

The focal plane array (FPA) on SAFARI'S band 3 (110/im- 
210/rai) is a 20 x 20 matrix, with a 1.9' x 1.9' field of 
view, so to map the entire 1.63° x 1.63° area we require 
(1.63 x 60/1. 9) 2 ~ 2600 pointings of the instrument. 

The 5— er, 1-hour sensitivity of SAFARI at 158 /xm is 
So =1.6 10" 19 W.m" 2 . Therefore, the time t(I) needed to 
detect an emerging fluxl is tq(Sq /I) 2 , with tq = 5 hours. 
This allows us to estimate the time needed for each of the 
2600 pointings, which ranges from ~1 to ^24 seconds. 
By adding these up, the total time needed to map the 
entire sky area is T ~ 4.5 hours. For the same sky area, 
Herschel would require some 900 hours, making such a 
project unrealistic. 

Of course this neglects the time needed to scan the 
band, since in Fourier- Transform Spectrometers such as 



SAFARI, this is done by translating a mirror. Depending 
on the spectral resolution R required, this can take from 
~ 10 s (R = 100) to up to 3.5 min (R = 2000). For the 
2600 pointings mentioned above, this means respectively 
7.2 hours and 152 hours. These are long overhead times, 
but one gets the full 34-210 /im band in return, which may 
contain many other interesting lines, such as [Nil] at 122 
/im and 206 /im. 

6. Conclusions and perspectives 

This work demonstrates the ability of SPICA/SAFARI to 
map large areas of the sky in the [Cll] line in a much 
shorter time than what would be possible with Herschel, 
making it clear that the extraordinary sensitivity of the 
proposed mission is a definitive asset regarding this type 
of project. 

On another level, the work presented here is part of an 
ongoing ASTRONET project dubbed STAR FORMAT, 
which is a German-French collaboration whose two main 
objectives are : 

1 / Producing databases to publicize simulation results 
from MHD and PDR codes; 

2/ Developing MHD, PDR and radiative transfer codes 
in an interoperable fashion, to allow for a much improved 
treatment of the physics of the ISM, especially regarding 
the formation of molecular clouds and dense cores. 

With respect to this project, our work shows that a 
truly three-dimensional PDR code is a must, to adequately 
deal with the complex geometries of the structures formed 
in the MHD simulations. The development of parallel com- 
puting schemes is also quite inevitable, and it should be 
noted that the PDR code can already be run on a grid. 
It is hoped that the present paper can be used as ground 
work for these developments. 
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